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Abstract: The shear viscosity coefficient of strongly coupled boundary gauge theory plasma 
depends on the horizon value of the effective coupling of transverse graviton moving in black 
hole background. The proof for the above statement is based on the canonical form of gravi- 
ton's action. But in presence of generic higher derivative terms in the bulk Lagrangian the 
action is no longer canonical. We give a procedure to find an effective action for graviton 
(to first order in coefficient of higher derivative term) in canonical form in presence of any 
arbitrary higher derivative terms in the bulk. From that effective action we find the effective 
coupling constant for transverse graviton which in general depends on the radial coordinate 
r. We also argue that horizon value of this effective coupling is related to the shear viscosity 
coefficient of the boundary fluid in higher derivative gravity. We explicitly check this pro- 
cedure for two specific examples: (1) four derivative action and (2) eight derivative action 
(Weyl 4 term). For both cases we show that our results for shear viscosity coefficient (upto 
first order in coefficient of higher derivative term) completely agree with the existing results 
in the literature. 
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1. Introduction 

The AdS/CFT correspondence is a powerful tool to study different properties of strongly 
coupled gauge theory in terms of dual (super) gravity theory in AdS space. In low frequency 
limit the boundary field theory can be described by hydrodynamics. In this limit different 
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transport coefficients like shear viscosity, diffusion constant, thermal and electrical conduc- 
tivity of strongly coupled boundary fluid have been computed in the context of AdS/CFT 
(see [1] - [31]). 

In [1], the authors evaluated the shear viscosity coefficient of boundary fluid using Kubo 
formula. This formula relates the shear viscosity to two point function of energy momentum 
tensor in zero frequency limit. On the other hand from field operator correspondence of 
the AdS/CFT conjecture we know that energy momentum tensor of boundary field theory 
is sourced by bulk graviton excitations. Therefore in the context of AdS/CFT, to calculate 
thermal two point correlation function of field theory energy momentum tensor we need to 
add small perturbations to the bulk metric. In [1], the authors considered graviton excita- 
tions polarized parallel to the black brane (i.e., xy components are turned on) and moving 
transverse to it. When one sends the gravitons to the brane, there is a probability that it will 
be absorbed by the brane. They calculated the absorption coefficient and showed that it is 
related to two point functions of energy momentum tensor of boundary fluid. 

To calculate the absorption coefficients, one needs to solve the wave equation for trans- 
verse gravitons. In presence of any higher derivative terms in the bulk action the solution may 
be technically difficult in general [32, 33, 34]. Recently there is a proposal that the shear vis- 
cosity of strongly coupled boundary gauge theory plasma is related to the effective coupling of 
graviton calculated at the black hole horizon [35, 36]. In [37], using membrane paradigm, the 
authors have confirmed that at the level of linear response the low frequency limit of strongly 
coupled boundary field theory at finite temperature is determined by the horizon geometry of 
its gravity dual. They have proved that generic boundary theory transport coefficients can be 
expressed in terms of geometric quantities evaluated at the horizon 1 . In particular, they have 
found that the shear viscosity coefficient is given by transverse graviton coupling computed 
at the horizon. The novelty of this result is that one does not need to solve the equation of 
motion for the graviton to calculate the thermal Green function. From graviton's action one 
can easily read off the coupling constant and hence determine the shear viscosity coefficient. 

To find the effective coupling of gravitons one has to find the general action. This can be 
achieved in the following way. Consider the Einstein-Hilbert action with negative cosmological 
constant 

/ = T6k/ dW ^ (ii + 12) - <L1) 

The equation of motion obtained from this action has a black hole solution. We denote this 
background solution by g$ . Now we consider fluctuation about this spacetime in xy (for 
example) direction 2 , 

9xy = g { $ + e h xy (r, x) = gVj (1 + e x)) . (1.2) 
Then substituting the metric with fluctuation in the action (1.1) and keeping terms up to 
^ee [38] also. 

2 Notations: x denotes the boundary coordinates, x = {t, x). 
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order e 2 



we get the action for graviton. The form of this action is, 



S ~ 



16vrG 5 



1 



/ 



d 4 k 



dr (a{r)<p'{r, k)(f)'(r, —k) + b(r)c/)(r, k)<p(r. 



k)) 



(1.3) 



where, 




k = {—uj,k} and ' ' ' denotes derivative with respect to r. The effective coupling is related 
to the coefficient of 4> 2 i.e., a (we have reviewed this calculation in section 2). 



This gives the correct viscosity coefficient for the Einstein-Hilbert gravity. But it is not 
obvious how to generalize this approach for higher derivative case. The proof given in [37] 
was based on the canonical form (1.3) of graviton's action. In presence of arbitrary higher 
derivative terms in the bulk, the general action for the perturbation h xy does not have the 
above form (1.3). Rather it will have more than two derivative (with respect to r) terms. 
[37, 39] have considered Gauss-Bonnet term in the bulk action. In general, presence of R a bR ab 
and R a i, C( iR abcd terms in the bulk result terms like <p" 2 and (j>'</>" in the action for h xy . For 
Gauss-Bonnet combinations these terms get canceled and the general action still has the form 



In this paper we have considered generic higher derivatives terms in the bulk Lagrangian. 
We have given a procedure to construct an effective action S cS for transverse graviton of 
the form (1.3) in presence of any higher derivative terms in the bulk. The details of the 
construction is given in section (3). Our construction ensures that in low frequency limit, 
the calculations of retarded Green function (imaginary part) using either effective action or 
original action are same. Therefore following the similar argument given in [37], we can 
relate the shear viscosity coefficient of the boundary fluid with the horizon value of the 
effective coupling obtained from S eS (section 4). In section (5) we have also discussed how 
membrane fluid captures the properties of boundary fluid in low frequency limit in generic 
higher derivative gravity. We have checked our procedure for two cases: 

• General four derivative terms, (section (6)) 

• Weyl A term which arises in type II string theory (section (7)). 

In both examples we get exact agreement between our results and the results that already 
exist in the literature [33, 34, 40]. Hence we conclude that: 

The shear viscosity coefficient of the boundary fluid is given by the horizon value of the effec- 
tive coupling of transverse graviton obtained from its effective action in presence of arbitrary 
higher derivative terms in the bulk. 



(1.3). 
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2. Shear Viscosity from Effective Coupling 



In this section we briefly review how to calculate the shear viscosity coefficient of the boundary 
fluid from the effective coupling constant of transverse graviton in Einstein-Hilbert gravity 

We first fix the background spacetime. We start with the following Einstein-Hilbert 
action in AdS space. 

/= — — ! d b x^g{R + 12) . (2.1) 



16ttG 5 _ 

Here we have taken the radius of the AdS space 1. The background spacetime is given by the 
following metric 3 

dr 2 1 

ds 2 = -h t (r)dt 2 + -f— + -dx 2 (2.2) 
h r (r) r 



where, 



and 



1 -r 2 

h t( r ) = ■ ( 2 - 3 ) 



h r (r) =4r 2 (l-r 2 ) . (2.4) 
The black hole has horizon at tq = 1 and the temperature of this black hole is given by, 

T=-. (2.5) 

7T 

We consider the following metric perturbation, 

9x« = 9 { S + Ky{r,x) = gWj(l + e$(r,x)) (2.6) 

where e is an order counting parameter. We consider terms up to order e 2 in the action of 
<J?(r, x). The action (in momentum space) is given by, 



1 f du;d 3 k 

-dr 



Ai,i(r)<i/(r,-k)<l/{r,k) 

(2.7) 



16vrG 5 J (2vr) 4 

+Ai, (r, k)<f>(r, -k)<j)'(r, k) + Ao. (r, k)<f>{r, k)<f>(r, -k) 



where, Aij(r, k) are functions of r and k and (f>(r,k) is given by (1.4). Up to some total 
derivative the action (2.7) can be written as 4 , 



°^dr (A?\r)4>'(r, -k)<f>'(r, k) + Af(r, fc)0(r, k)<f>(r, -k)) (2.8) 



dujd 3 k 
16ttG 5 J ~(2i 

where, 

Ai>{r) = - (2.9) 



,(0)^ r z -\ 



3 We are working in a coordinate frame where asymptotic boundary is at r — ► 0. 

4 Though throughout this paper we have written the four vector k, but in practice we have worked in k — > 
limit. In all the expressions we have dropped the terms proportional to k or its power. 
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and 

4 0) M) = 4r2( f_ r2) • (2-10) 

This can be viewed as an action for minimally coupled scalar field 4>(r, k) with effective 
coupling given by, 

K -M = i«S7=^- (2 ' u) 

Therefore according to [37, 39] the effective coupling K e R calculated at the horizon tq gives 
the shear viscosity coefficient of boundary fluid, 

r l = r~ l (-2K eS (r Q )) 
= — !— . (2.12) 

3. The Effective Action 

Having understood the above procedure to determine the shear viscosity coefficient from the 
effective coupling of transverse graviton it is tempting to generalize this method for any higher 
derivative gravity. As we discussed in the introduction, the first problem one faces is that the 
action for transverse graviton no more has the canonical form (2.7). For generic 'n' derivative 
gravity theory the action can have terms with (and up to) 'n' derivatives of <J>(r, x). Therefore, 
from that action it is not very clear how to determine the effective coupling. In this section 
we try to address this issue. 

We construct an effective action which is of form (2.8) with different coefficients capturing 
higher derivative effects. We determine these two coefficients by claiming that the equation 
of motion for <p(r, k) coming from these two actions (general action and effective action) are 
same up to first order in perturbation expansion (in coefficient of higher derivative term). 
Once we determine the effective action for transverse graviton in canonical form then we 
can extract the effective coupling from the coefficient of (f)'(r, k)(f)'(r, —k) term in the action. 
Needless to say, our method is perturbatively correct. 

3.1 The General Action and Equation of Motion 

Let us start with a generic 'n' derivative term in the action with coefficient fi. We study this 
system perturbatively and all our expressions are valid up to order n. The action is given by, 

where, IZ^ is any n derivative Lagrangian. The metric in general is given by (assuming 
planar symmetry), 

ds 2 = -(h t (r) + fi hl n \r))dt 2 + - ^ + V h^(r))dx 2 (3.2) 

h r (r) + fi hi '(r) ' 
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where h[ n \ h!^ and are higher derivative corrections to the metric. 

Substituting the background metric with fluctuations in the action (3.1) (we call it general 
action or original action) for the scalar field cf)(r, k) we get, 

1 r rl^h n 

S = lfoFGe J ^ k) ^ P){r > - k ^ (q) ( r > k ) (3-3) 

where, <f>^ (r, k) denotes the p th derivative of the field 4>(r, k) with respect to r and p + q < n. 
The coefficients A Pt q(r, k) in general depends on the coupling constant jjl. A pa with p + q > 3 
are proportional to fi and vanishes in fi — > limit , since the terms <f>( q ) with p + q > 3 
appears as an effect of higher derivative terms in the action (3.1). Up to some total derivative 
terms, the general action (3.3) can also be written as, 



1 f rrk 

S = -HTrr / ^rfrj^ir.fc)^^-*)^^*), n even 
167rG 5 J (2vr) 4 ^ 

= m^jj^Y dr S ^MW^-AO^M), n odd. (3.4) 



p=0 

The equation of motion for the scalar field (p(r, k) is given by, 



; () , ^(P)(r,fe) 

( "^ /2 / (i\ p a£({0( m )}) 

> — r I — ttt — \ = 0. n odd 3 - 5 ) 
^ V dry d(p(p)(r,k) 



n even 



p=0 

where £({c^ m )}) is given by 

£({0M}) = ^i p (r,fc)^(r,-fc)^(r,fc) . (3.6) 
p 

We analyze the general action for the scalar field 4>(r, k) and their equation of motion 
perturbatively and write an effective action for the field (f)(r, k). 

The generic form of the equation of motion (varying the general action) upto order fi is 
given by, 

Mr, k)<P(r, k) - Mr, k)<f>'(r, k) - A^r, k)^"(r, k) = ji £({<f> ip) }) + 0(p 2 ) (3.7) 

where f{{<t> {p) }) is some linear function of double and higher derivatives of (p(r,k), coming 
from two or higher derivative terms in action (3.3). The zero th order (// — > 0) equation of 
motion is given by, 

A { o\r, k)</>(r, k) - Al°\r, k)cp'{r, k) - Af\r, k)</>"(r, k) = (3.8) 
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where, A^ is the value of A p at n — > 0. From this equation we can write 4>"(r, k) in terms of 
4>'(r, k) and (f)(r, k) in \x — ► limit. 

0"(r, fc) = k ) ~ A m^ k) J ^ fc ) • ( 3 - 9 ) 

Af\r,k) A { i\r,k) 
Then the full equation of motion can be written in the following way, 

4 0) (r, k)cP(r, k) - X (0) (r, k)<f>'(r, k) - A? (r, fc)0"(r, fc) = /x jfy(r, fe), </»'(r, k), <j>" (r, k), ...) 

+0(^ 2 ) . (3.10) 

Since the right hand side of equation (3.10) is proportional to fx, we can replace the <f)"(r, k) and 
other higher (greater than 2) derivatives of (f)(r, k) by its leading order value (3.9). Therefore 
up to order fi the equation of motion for 4> is given by, 

A { o\r, k)(j)(r, k) - A} 0) (r, k)<f>'(r, k) - Af\r, k)<f>"{r, k) = n T{(f>(r, k), <f)'(r, k)) 

+eV) 

= (r,k)+F <f>(r,k)) 

+CV) (3.11) 

where JF and T\ are some function of r. This is the perturbative equation of motion for the 
scalar field <fi(r, k) obtained from the general action (3.3). 

3.2 Strategy to Find The Effective Action 

In this subsection we describe the strategy to write an effective action for the field </>(r, k) 
which has form (2.8) with different functions. The prescription is following. 

• (a) We demand the equation of motion for (f>(r, k) obtained from the original action 
and the effective action are same upto order fj,. This will fix the coefficients of 2 and 
(p 2 terms in effective action. 

Let us start with the following form of the effective action. 



dud?k 



,(o) 



{A ( °\r, k) + /i£i(r, k))4>'{r, -k)<j>'(r, k) 



+(Ar (r, k) + fiB (r, k))cj>(r, k)cf>(r, -k) 



(3.12) 



The functions B$ and B\ are yet to be determined. We determine these functions by 
claiming that the equation of motion for the scalar field <f)(r, k) obtained from this 
effective action is same as (3.11) up to order fj,. The equation of motion for (f>(r, k) from 
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the effective action is given by, 

A ( o\r, k)(f)(r, k) - A}° ] (r, k)<j>'{r, k) - Af ] (r, k)<f/'{r, k) 



= „ U(r, k) - A \^ k) B^ k)) J>>(r, k) 
V A?\r,k) J 



+f i (b 1 {t, k) A f n . {r ' k) - Bo(r, k)) <f>(r, k) + 0(fi 2 ) . (3.13) 



A?\r,k) 



Therefore comparing with (3.11) we get, 



K(r, fc) " )J ' > Bi(r, k) - ^i(r, fc) = (3.14) 



4 0) (r,A;) 
and 

B (r, fc) = Bi(r, *0 °/ ' J - Mr, k) . (3.15) 
A\ (r, fe) 

The solutions are given by, 

B 1 (r,k) =Af\r,k) [ dr + KA?\r,k) 

J A\ (r, A;) 

= Bi(r,fc) + ^S 0) (r,fc) (3.16) 

and 

B = B (r,k) + KA ( o ) (3.17) 
for some constant k. We need to fix this constant. 

• (b) Condition (a) can not fix the overall normalization factor of the effective action. In 
particular we can multiply it by (1 + fiT) (for some constant T) and still get the same 
equation of motion 5 . Considering this normalization, the effective action is given by, 



1 + n T f dwd 3 k 



(A?\r, k) + nB\ (r, k))<f>'(r, -fc)^'(r, k) 

(3.18) 



/" dud 

J w 

+ (4 0) (r, fe) + fJ,Bo(r, k))cf>(r, fc)0(r, -fc) 
Substituting the values of £>'s (3.16) and (3.17) we get, 

S ca = (1 + Mr + «))S (0) + // y ^ (Bi(r, fc)^'(r, -fc)^'(r, fc) + B (r, fc)0(r, -fc)0(r, fc)) 

(3.19) 

where is the effective action at fi — > limit. This implies that the integration 
constant k can be absorbed in the overall normalization constant T. Henceforth we will 
denote this combination as T. 

Our prescription is to take T to be zero from the following observation. 



'We are thankful to Ashoke Sen for raising this point. 
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The shear viscosity coefficient of boundary fluid is related to the imaginary part 
of retarded Green function in low frequency limit. The retarded Green function 
Gxy,xy(k) is defined in the following way. The on-shell action for graviton can be 
written as a surface term, 

/d 4 k 
j^Mk)Gxy,xy(k,r)M-k) (3.20) 

where 4>o(k) is the boundary value of cj)(r, k) and G xyxy is given by, 

G«y, xy (k) = Km2g xytXy (k,r) (3.21) 

and shear viscosity coefficient is given by 6 , 



77 = lim 



LLL ^xy,xy 



(k) 



(computed on — shell) . (3.22) 



— Now it turns out that the imaginary part of this retarded Green function obtained 
from the original action and effective action are same upto the normalization con- 
stant T in presence of generic higher derivative terms in the bulk action. Therefore 
it is quite natural to take T to be zero as it ensures that starting from the effective 
action also one can get same shear viscosity using Kubo machinery. To show that 
the above statement is true we do not need to know the full solution for </>, in other 
words to find the difference between the two Green functions one does not need to 
calculate the Green functions explicitly. Assuming the following general form of 
solution for (f) 

0~ (1 -r 2 )~ iul(i (l + iuj^ui(r)) (3.23) 
it can be shown generically. In appendix A we have given the proof. 

— Because of the canonical form of the effective action, it follows from the argument 
in [37] and the statement above, that the shear viscosity coefficient of boundary 
fluid is given by the horizon value of the effective coupling obtained from the 
effective action in presence of any higher derivative terms in the bulk action. We 
discuss elaborately on this point in section (4). 

(c) After getting the effective action for </>(r, k), the effective coupling is given by, 

_ i (3 . 24) 



where g rr is the Vr' component of the inverse perturbed metric and \/—g is the deter- 
minant of the perturbed metric. Hence the shear viscosity coefficient is given by, 



3 



r, = r 2 (-2K eS (r = r )) (3.25) 
where rg is the corrected horizon radius. 



6 To calculate this number one has to know the exact solution, i.e., the form of £ and the value of (3 in 
(3.23). 
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To summaries, we have obtained a well defined procedure to find the correction (up to order 
fi) to the coefficient of shear viscosity of the boundary fluid in presence of general higher 
derivative terms in the action. 



4. Flow from Boundary to Horizon 



Following [37], let us define the following linear response function 

where II(r, fc) is conjugate momentum of the scalar field (f> (with respect to a foliation in the 
r direction), 

n(r, fc) = (Af\r, k) + fiBi (r, fc)) </>'(r, —k) 



(4.2) 



where K e s{r) = IQirG^K^r). Now we will show, using the equation of motion, that the 
function II(r, k) and the combination u<f>(r, k) is independent of the radial coordinate r in 
k — > limit. The equation of motion is given by, 

d ~ 



dr 



(A?\r, k) + fiB^r, fc)) 4>'(r, fc) = {j^\r, fc) + fiB (r, fc)) 0(r, fc) 



d_ 

dr 



n(r, fc) 



(4 O) (r,fc)+^ o (r,fc))0(r,fc) . (4.3) 



Since A^ ~ to 2 , therefore it follows from (4.3) and (4.2) that, in \x — > limit II (r, fc) and 
uj(j)(r, fc) are independent of r. But this is true even in jjl ^ case. To understand this we 
note that, function Aq in (3.4) is proportional to to 2 in general . Therefore it follows from 
(3.9), (3.11) and (3.15) that Bq is also proportional to uj 2 . Hence, in presence of higher 
derivative terms also it follows from (4.2) and (4.3) that the function II(r, fc) and w(f)(r, fc) are 
independent of radial direction r in low frequency limit. 

Therefore this response function calculated at the asymptotic boundary and at the hori- 
zon gives the same result and is equal to the shear viscosity coefficient. One can calculate 
the function \ an d it turns out that, 

IIIU" 

X(r = 0, fc 



()s = ImG xy,o ■!, 



IUJ 



X(r = r ,k -> 0) 



r Q 3/2 A?\r,k) + fiB!(r,k) 



r Q 2 (-2K cS (r )) 



(4.4) 



8vrG 5 V=99 rr 

Thus, shear viscosity coefficient of boundary fluid is related to horizon value of graviton's 
effective coupling obtained from the effective action. 

7 In general when we write action (3.4) action (3.3) we get some terms like u> 2 (j) 2 + Z(r)<f> 2 . The function 
Z(r) is zero when background equation of motion is satisfied. We have explicitly checked this for two, four 
and eight derivative case. 
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5. Membrane Fluid in Higher Derivative Gravity 

The UV/IR connection tells us that the boundary field theory physics in low frequency 
limit should be governed by the near horizon geometry of its gravity dual. In [37], the 
authors have established a connection between horizon membrane fluid and boundary fluid in 
linear response approximation. They considered a mass less scalar field (with action given in 
(2.8)) outside the horizon and studied the response of the membrane fluid to this bulk scalar 
field. They defined a membrane charge H m b which is equal to the conjugate momentum of 
the scalar field 4> (with respect to a foliation in the r direction) at the horizon. Imposing 
regularity condition on the scalar field at the horizon they interpreted the membrane charge 
Il m b as a response of the horizon fluid to the scalar field. Considering the scalar field cj> to be 
bulk graviton excitation (h%), Il m b gives the shear viscosity of the membrane (horizon) fluid 
which is also equal to horizon value of the effective coupling of graviton. In this way, they 
proved that the shear viscosity of boundary fluid is related to that of membrane fluid. 

In higher derivative gravity, since the canonical form of the action (2.8) breaks down, it 
is not very obvious how to define the membrane charge Il m b. Instead of the original action if 
we consider the effective action (3.12) for graviton then it is possible to write the membrane 
action perturbatively and define the membrane charge (II m b) in higher derivative gravity. As 
if the membrane fluid is sensitive to the effective action 5 efl in higher derivative gravity. 

Following [37] we can write the membrane action and charge in the following way (in 
momentum space) 

where a^ u is the induced metric on the membrane and II(r, k) is given by (4.2) and the 
membrane charge is given by, 

n mb = = -K^Jg^dr^r, k)\ rQ . (5.2) 

Imposing the in-falling wave boundary condition on 0, it can be shown that the membrane 
charge II m b is the response of the horizon fluid to the bulk graviton excitation and the 
membrane fluid transport coefficient is given by, 

Vmb = K cS (r ) . (5.3) 

Hence, we see that even in higher derivative gravity the shear viscosity coefficient of 
boundary fluid is captured by the membrane fluid. 

6. Four Derivative Lagrangian 

In this section we apply our effective action approach to calculate the correction to the shear 
viscosity in presence of general four derivative terms in the action. The four derivative bulk 
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action we consider is of the following form 

S = J^r^r J d 5 x [R + 12 + n [c x R 2 + c 2 R ab R ab + c 3 R abcd R abcd ^j 

with constant c±, c 2 and C3. The background metric is given by, 

1 



- f -^dt 2 + dr 



4r 2 f(r] 



+ -dx 2 



where, 

The position of the horizon is given by, 



f(r) = 1 - r 2 + |(4(5ci + c 2 ) + 2c 3 ) + 2^c 3 r 4 



/M = 



which implies that, 



ro = 1 + -(5ci + c 2 + 2c 3 )/* + O^ 2 ) 



The temperature of this black hole is given by, 

1 (5ci + c 2 - 7c 3 )/i 



T= - + 

7T 



3vr 



+ O {^) 



In this coordinate frame the boundary metric is given by, 

ds\ = (-f(0)dt 2 + dx 2 ) 



(6.1) 

(6.2) 
(6.3) 

(6.4) 
(6.5) 

(6.6) 

(6.7) 



which is not Minkowskian. Therefore we rescale our time coordinate to make the boundary 
metric Minkowskian. We replace, 

t 



V7W) 



in the metric (6.2). The rescaled metric is, 



ds 2 = - 



f(r) 



dt 2 + 



dr 2 



+ -dx 2 . 



f(0)r~" ' 4r 2 f(r) r 
This is our background metric and we consider fluctuation around this. 

6.1 The General Action 

In this theory, the general action for the scalar field 4>(r, k) is given by, 



(6.8) 



(6.9) 



1_ f d 4 k 



16vrG 5 



■dr 



Af B (r, k)(f)(r, k)<f>(r, -k) + Af B {r, k)<f/(r, k)<f/{r, -k) 
+A% B (r, k)4>"{r, k)<f>"(r, -k) + A^ B (r, k)<j>{r, k)4>'{r, -k) 
+A% B (r, k)<f>(r, k)<f>"(r, -k) + A% B (r, k)<f>'(r, k)<p"(r, —k) 



(6.10) 
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where the expressions for Af B s are given in appendix B. Up to some total derivative terms 
this action can be written as, 

1 f rl A k r 

S = J J^dr A^^r, k)<f>(r, —k) + Af B cj ) '(r, k)<//{r, -k) + A° B cf>"(r, k)<ft'(r, -k) 

(6.11) 

where, 



A$ B = A^(r,k) 



GB, 



A'f B (r,k) Al GB (r,k) 



+ 



A? B = A° B (r,k)-Ag B (r,k) 
Af B = Af B (r,k). 
6.2 The EfTective Action and Shear Viscosity 



A'f B (r,k) 



(6.12) 



Following the general discussion of section (3) we write the effective action for the scalar field, 



GB = (1 + I>) 
off 16vrG 5 



f d A k 

J WT 4 



(4 0) (r, k) + nB? B (r, k))<t>'(r, -k)4>'{r, k) 



+ (A^"(r, k) + ^Bq B (r, k))4>(r, k)<j>{r, -k) 



(6.13) 
(6.14) 



To evaluate the functions Bf B and Bq B and to fix the normalization constant T, we follow 
the strategy given in section (3.2). Comparing the equation of motion for <fi(r,k) from two 
actions we get the function Bf B and Bq B of the following form, 

,2 

(I0(llr 2 - 13)ci + (22r 2 - 26)c 2 + (11 - 25r 2 + 6r 4 )c 3 ) 



B GB 



12r 2( 1 _ r 2)2 



B? B = — ((110 - 130r 2 )ci + (22 - 26r 2 )c 2 - (13 - 23r 2 + 18r 4 )c 3 ) . (6.15) 

The normalization constant T = (appendix A). 

Now we can calculate the effective coupling using the formula (3.24). It turns out to be, 

K cS (r) = + ( 4 ( 5c i + C2 ) " 2 ( X " m) • (6-16) 

Therefore the shear viscosity is given by, 

^ = (- 2 ^ff(ro)) 



1 



1 



16vrG 5 r 3 / 2 



(1 - 8(5 Cl + c 2 )fi) 



1 



16ttG 5 

This result is in agreement with [33, 34, 41]. 



(1 - 9(i (5ci + c 2 ) - 2fi c 3 ) 



(6.17) 
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7. String Theory Correction to Shear Viscosity 



In this section we apply the effective action approach for eight derivative terms in the La- 
grangian. We consider the well known Weyl 4 term. This term appears in type II string 
theory. The five dimensional bulk action is given by, 



s -mkJ d ^( R + 12 + " wi ") 



where, 

Ti/(4) /shrank/- 

and the weyl tensors C a b c d are given by, 



tia(4) sihmnksi r 1 rs Pr<1 i r<hkmnr< C< rsp^q 

<■ milium /, ^ rsk 2 ^pqmn^h ' 



rsk 



Cabal — Rabcd + ^(dadRcb + 9b C R a d ~ 3acRdb ~ 9bdRca) + — (9ac9bd ~ 9ad9cb)R ■ 

The background metric is given by [42, 43], 

ds 2 = - K ' (1 + 45^r 6 - 75^r 4 - 75^r 2 ) dt 2 

+ 4 ^ 1 2 (l - 285^r 6 + 75,ur 4 + 75/ir 2 ) dr 2 + ^-dx 2 . 

The temperature of this black hole is given by, 

T = - (1 + 15/x) . 

7T 



(7.1) 



(7.2) 



(7.3) 



(7.4) 



(7.5) 



The horizon is located at tq = 1. 
7.1 The General Action 

Putting the perturbed metric in (7.1) we get the general action for the scalar field (p(r, k), 

d 4 k 



-dr 



Af{r, k)(f)(r, k)(p(r, -k) + Af {r, k)</>'(r, k)4>'{r, -k) 



16ttG 5 J (2vr) 4 

+Af(r, k)(j)"{r, k)^"{r, -k) + Af{r, k)<f>(r, k)4>'{r, -k) 

+Af{r, k)<j)(r, k)(f)"(r, -k) + Af (r, k)<f>'(r, k)<f/'(r, -k) 



(7.6) 



The coefficients A^s are given in appendix (C). Like four derivative case, up to some total 
derivative terms this action can be written as, 



1 



d A k 



16vrG 5 J (2vr) 4 



dr 



Atf&ir, k)<f>(r, -k) + Af c/>'{r, k)</)'(r, -k) + A?<i/'(r, k)</)"(r, -k) 



(7.7) 
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where, 



w i 



A'W(r,k) Al w {r,k) 



+ 



= A™(r,k)-A™(r,k)- 



. w 



W i 



A' 6 w (r,k) 



A? =A?(r,k). 
7.2 The Effective Action and Shear Viscosity 

We write the effective action for the scalar field in the following way, 
(l + r» f d A k 



S. 



w 

eff 



16vrG 5 



f d*k 

J W) 



(Af\r, k) + nB^ir, k))<f/(r, -k)<f>'(r, k) 



on 

+(4 0) (r, k) + ^(r, k))<f>(r, k)<j>(r, -k) 
The functions BY and BY are given by, 



BY(r,k) 



(663r 6 - 573r 4 + 75r 2 ) 



BY(r,k) 



4r 2 (r 2 - 1) 
(r 2 - l) (I29r 6 + 141r 4 - 75r 2 ) 



The normalization constant T = (Appendix A). 

The effective coupling constant is given by (3.24), 

1 A { i ] (r,k)+nBY(r,k) 



K cS (r) 



16vrG 5 



99' 



16ttG 5 

Therefore the shear viscosity is given by, 



L- (-i (1 + 36,^(6-^))) 



16ttG 5 

and shear viscosity to entropy density ratio 



r? = r >(-2K cS (r )) 

1 (1 + 180 fj.) , (r = 1) 



5 = J- (1 + 120 m) 

S 47T 



where entropy density s is given by [42, 43], 

1 



s = 



4G 5 



(1 + 60 a*) • 



These results agree with [40] 8 . 



(78) 



(79) 
(7.10) 

(7.11) 
(7.12) 



(7.13) 



(7.14) 



(7.15) 



(7.16) 



3 In fact, in [32] the result for rj/s was not correct. Later the author(s) corrected their results in [40]. 
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8. Discussion 



We have found a procedure to construct an effective action for transverse graviton in canonical 
form in presence of any higher derivative terms in bulk and showed that the horizon value of 
the effective coupling obtained from the effective action gives the shear viscosity coefficient 
of boundary fluid. Our results are valid upto first order in fi. We discussed two non trivial 
examples to check the method. We have considered four derivative and eight derivative 
{Weyl 4 ) La grangian and calculated the correction to the shear viscosity using our method. We 
found complete agreement between our result and the results obtained using other methods. 

Since the equation of motion for scalar field 4>(r, k) obtained from effective and original 
actions are same, these two actions should be related by some field re-definition. If one finds 
such field re-definition then the normalization of the effective action will be fixed automati- 
cally 9 . 

In [35] the authors have proposed a formula for shear viscosity for generalized higher 
derivative gravity in terms of some geometric quantity evaluated at the event horizon (like 
Wald's formula for entropy). Though their proposal gives correct results for Einstein-Hilbert 
and Gauss-Bonnet action but unfortunately we are unable to get the correct result for Weyl 4 
term. We find the shear viscosity coefficient for Weyl 4 term (using their proposal) 

"=16k (1 + 20 " ) 

which implies, 

2 = i- (1 - 40^) . (8.2) 

These issues are under investigation [44]. 

In this paper we have concentrated on a particular transport coefficient, namely the 
shear viscosity coefficient. But the other transport coefficients like electrical and thermal 
conductivity of boundary fluid can also be captured in terms of membrane fluid. It would 
also be interesting to study these other transport coefficients in the context of higher derivative 
gravity. 
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Appendix 



A. Fixing the Normalization Constant 



In this appendix we fix the normalization constant T. We consider a general class of action 
for <fr which appears when the higher derivative terms are made of different contraction of 
Ricci tensor, Riemann tensor, Weyl tensor, Ricci scalar etc. or their different powers. Since, 
all these tensors involve two derivatives of metric they can only have terms like d a d^{r,x) 
and its lower derivatives. Therefor the most generic quadratic (in 3>(r, x), in linear response 
theory) action for this kind of higher derivative gravity has the following form (in momentum 
space) 10 



a _ i r d*k A 



al(r)(/>(r) 2 + a2(r)c//(r) 2 + a4(r)<£(r)0'(r) 



+ji &6(r)4>" (r)4>' (r) + ji &3(r)4>"(r) + &5(r)(/)(r)(/)" (r 



(A.l) 



where, 



-8r 2 + uj 2 r + 



al(r 
a2(r 
a4(r 
a5(r 

and a3(r), a6(r), j'2(r), g2(r), h2(r 
Now let us write the effective Lagrangian as follows, 



4^3 ^5 

3 

= -3r H h H h2(r) 

r 

6 

= -^-2 + Mg2(r) 
4 

= -4r + - + (i j2(r) 



+ /x f2(r) 



(A.2) 



and f2(r) depends on higher derivative terms in the action. 



_ i + ^r y_d*fc_ 



16ttG 5 



>'• r 4r ( r 2 _ X y 0/( r )2 _ w 2^( r )2 

)4™ y 4r 2 (r 2 - 1) 

'/„\2 



+/z b2(r)0(r)^+ bl(r)0'(r) 



(A.3) 



From condition (a) of section (3) the solutions for 61 and 62 are given by, 

1 



61(r) 



; ((-4r 3 - 12r + w 2 )a3(r 



2r (r 2 - l) 2 

+ (r 2 - 1) (2^r 4 - a6'(r)r 3 - 4/tr 2 + 2a3'(r)r 2 

+2 (r 2 - 1) h2(r)r - 2 (r 2 - l) ]2{r)r + a6'(r)r + 2k + 2a3'(r))) 



(A.4) 



■'In all the expressions we have omitted k dependence of . 
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62(r) = - 



-{{uo A + U4r 3 uj 2 ) a3(r 



16r 2 (r 2 - l) 4 

+4 (r 2 - 1) (-4r 2 f2(r) (r 2 - l) 3 + (2r 2 g2'(r) (r 2 - l) 2 
+ (c^ 2 k - 2r 2 (r 2 - 1) j2») (r 2 - l) 
+rcj 2 a3"(r)) (r 2 - l) + (l - llr 2 ) K; 2 a3'(r))) . 



(A.5) 



The boundary term coming from the original action (after adding Gibbons-Hawking boundary 
terms) are given by, 



d A k 



16ttG 5 J (2vr 



+ 4>(r) 2 + r<f>'(r)(/>(r) - 



0'(r)0(r) 



+/U Qg2(r)0(r) 2 -ij2'(r)0(r) 2 

+h2(r)<^'(r)^(r) - j2(r)0'(r)^(r) - -a6' (r)^'(r)^(r) 
a3'(r) (0(r)w 2 + 4 (r 4 - l) <jj(r)) <f>(r) 

4 r ( r 2_ 1 ) 2 

a3(r) (6r0(r)w 2 + (r 2 - l) (8r 3 + 24r - w 2 ) <f/(r)) <f>(r) 

4r (r 2 — l) 3 
a3(r)0'(r) (0{r)uj 2 + 4 (r 4 - l) <j>'{r)) 



+a3(r)^(r) 



4r (r 2 — l) 2 

0(r)u 2 (r A -l)cp'{r) 



\ 2r(r 2 -l) 2 r(r 2 -l) 2 , 
And the boundary terms coming from the effective action are given by, 
1 f d 4 k 



(A.6) 



off 16itG 5 J (2vr) 4 
+ 



r \ 4>{r)<p'{r) 



2r(r 2 - I)' 



^ (f)(r)(2T (r 2 - l) 3 + (-a6'(r> 



+2a3'(r)r 2 + 2 (r 2 - l) h2(r)r - 2 (r 2 - l) j2(r)r 

+a6'(r)r + 2a3'(r)) (r 2 - l) + (-4r 3 - 12r + uj 2 ) &3(r))<t>' (r) 



(A.7) 



Let the form of the solution of 4> is given by, 

(l-r 2 )^ (l + i/3unF(r)) 



with 



F(0) = 0. 



(A.8) 



(A.9) 
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The imaginary part of the retarded Green function for original action is given by, 



-Im[G£ ,xy J 



lim 



[HP(4 (r 2 + 3) a3(r)r 2 + (r 2 



-2/3 + T , , 

r (r 2 — 1) \ 

-l)(F'(r)r 6 + a6'(r)r 4 - 3F'(r)r 4 - 2a3'(r)r 3 - 2 (r 2 - l) h2(r)r 2 
+2 (r 2 - 1) j2(r)r 2 - a6'(r)r 2 + 3F'(r)r 2 - 2a3' (r)r - F'(r))) 

(A.10) 



and imaginary part of the retarded Green function for effective action is given by 



lwr , - R(effective) i 

un l^- r xy,xy \ 



UJ 



lim 

r->0 



2/3 . 

V(r 2 -1) ; 

(2rr 4 - a6'(r)r 3 - 4rr 2 + 2a3'(r)r 2 + 2 (r 2 - l) h2(r)r 

-2 (r 2 - l) j2(r)r + a6'(r)r + 2r + 2a3'(r)) - 4r (r 2 + 3) a3(r 



-rF'(r) + 



F'(r) 







(A.ll) 



Therefore, in low frequency limit the difference between the imaginary part of retarded Green 
function coming from this two boundary terms are given by , 



lim — Im 

co^O UJ 



G 



jj(origmal) 

xy,xy 



— Im 

UJ 



G 



^(effective) 

xy,xy 



2n r . 



(A.12) 



Therefore for this general class of theory, 



r = . 



(A.13) 



The other kind of higher derivative theory one can consider is covariant derivatives acting 
on curvature tensors. In that case one can have a more general action like (3.3). For this 
kind of action the boundary terms one gets are of the form <p^(j)^ (here (p^ means n-th 
derivative of <j) with respect to r). Using the form of (f> given in (3.23) it can be shown that 
except (jy^cj) kind of terms, other boundary terms do not contribute in low frequency limit. 
For example, if we consider C n (j)^ 2 term in the original action, the the relevant boundary 
term which will contribute in low frequency limit is (— l)( ra+1 )(C n (/>( n ))( n-1 )0. One can check 
that though we need to add Gibbons- Hawking terms to make the variation of the action well 
defined but most of them are zero in low frequency limit. We have checked it for few nontrivial 
terms like, (p^ 2 and (p^ 2 and T turns out to be zero. But we expect it is true in general. 
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B. Expressions for A GB 



A ? B ^ k) = s i;^ r J[) - (l2FF^TF)) ((10ci(88r4 " n " V " 176r2 + 13uj2r + 88) 

+c 3 (144r 8 - 288r 6 + 66u; 2 r 5 + 232r 4 + 25wV - 4(3w 4 + 44)r 2 + 13u; 2 r + 88) 
+c 2 (176r 4 - 22w 2 r 3 - (3a; 4 + 352)r 2 + 26w 2 r + 176))//) + 0(// 2 ) 

r 

(10ci(13r 2 - 11) + 2c 2 (2r 4 + 17r 2 - 9) + c 3 (34r 4 + 9r 2 - 8cj 2 r - 3))/x 2 , 
H - h U(n ) 

Af B (r,k) = 4(c 2 + 4c 3 )r(r 2 -l) 2 // + 0(// 2 ) 

+ ^tA r(2(10ci(13r 4 + 20r 2 - 33) + c 2 (26r 4 + 3wV + 40r 2 + 3cA - 66) 

3r z (r A — 1) 

+c 3 (90r 6 - 89r 4 + 30u; 2 r 3 + 32r 2 + 6w 2 r - 33))//) + O (// 2 ) 

A as (r , k) _ _^21 

2(20ci(13r 2 - 11) + 2c 3 (18r 4 + r 2 - 11) + c 2 (52r 2 + 3c^ 2 r - 44))// , 2 \ 
+ 3r + U W ) 

A% B (r, k) = 8 (r 2 - l) (c 2 r 2 + 4c 3 r 2 + c 2 ) // + O (// 2 ) . (B.l) 
C. Expressions for A w 



A W = 8r 2 - (V 2 r 



4 r 3 (j.2 _ lj 

(-360r 9 - 240r 7 + 129w 2 r 6 + 1560r 5 - 300aA 4 + 8 (to 4 - 120) r 3 + 75lo 2 ) 7 , 9N 

+" w^w " " ~ + W 

3 ( r 2 — 1) 

Af = y — '- + r (-419r 6 + 668r 4 - 24wV + 8r 2 - 225) // + O (// 2 ) 

A W = 32r 5 p-lJ^ + O^ 2 ) 

w _ 2 (r 2 + 3) 2 (2045r 8 - 4185r 6 - 26w 2 r 5 + 2140r 4 - 2w 2 r 3 + 75r 2 - 75) /z 2 n 
^4 = ^2 ^T] h ° ) 

4 fr 2 - l) 

Af = i '- - 4 (r (I45r 6 - 220r 4 + 2wV + 75)) // + O (// 2 ) 

A™ = 32r 4 (2r 4 - 3r 2 + l) // + O (// 2 ) . (C.l) 
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